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Abstract 

This paper is an overview of our works which are related to investigations of the 
integrability of natural Hamiltonian systems with homogeneous potentials and New- 
ton's equations with homogeneous velocity independent forces. The two types of in- 
tegrability obstructions for these systems are presented. The first, local ones, are re- 
lated to the analysis of the differential Galois group of variational equations along a 
non-equilibrium particular solution. The second, global ones, are obtained from the 
simultaneous analysis of variational equations related to all particular solutions belong- 
ing to a certain class. The marriage of these two types of the integrability obstructions 
enables to realise the classification programme of all integrable homogeneous systems. 
The main steps of the integrability analysis for systems with two and more degrees of 
freedom as well as new integrable systems are shown. 

ke}rwords:integrability; differential Galois theory; Hamiltonian systems. 



1 Introduction 

It is hard to imagine yourself a world without integrable models. Teaching in such a world 
would be rather frustrating. If a theory has no solvable examples, it is difficult to explain 
that it is useful. Fortunately, in ours we have the harmonic oscillator — the guinea-pig which 
serves as pedagogical example for all theories. 

Integrable models are exceptional, but we do not neglect them. Still, as it was one 
century ago, finding a new non-trivially integrable model is a discovery. 

How we can find integrable systems? One way, which seems to be the most natural, is 
just a search them in the nature. That is, start with more or less general model, and look for 
some special cases when the model is integrable. The other way is to construct integrable 
systems. It appears that the first approach is much more difficult than the second one. The 
reason of this is obvious: we know only few general methods which give strongly enough, 
and computable necessary conditions for the integrability. 

In this paper we consider only classical dynamical systems which are described by or- 
dinary differential equations. There is no a unique definition of such systems and there is 
no a unique method for study of their integrability. Nevertheless, in this paper we will con- 
centrate mainly on applications of only one quite new theory. It was developed by Baider, 
Churchill, Morales, Ramis, Rod, Simo and Singer in the end of the XX century, see fill; Hi 



and references therein. In the context of Hamiltonian systems it is called the Morales-Ramis 



theory. It arose from very long searching for relations between the branching of solutions 
of differential equations considered as functions of complex time, and the integrability. In 
the context of Hamiltonian systems these relations were found by S. L. Ziglin [50; 51]. His 
elegant theory expresses necessary conditions for the integrability in terms of properties of 
the monodromy group of variational equations along a particular solution. The Morales- 
Ramis theory, in some sense, is an algebraic version of the Ziglin theory. It formulates the 
necessary conditions for the integrability in terms of the differential Galois group of the 
variational equations. 

During one and half decade the Morales-Ramis theory was applied successfully for 
study the integrability of numerous systems, see an overview paper I©]. Let us mention 
two biggest successes of this theory. It was applied to prove the non-integrability of the 
planar three body problem jlj; |47|: i49|] , and to prove the non-integrability of the Hill 
lunar problem [371]. Moreover, as it is well known, the first proof of the fact that the problem 
of the heavy top is integrable only in the classical cases was done by S. L. Ziglin in i5ll] 
and it is based on his theory. An alternative proof based on differential Galois approach is 
given in 

The above examples show the real power of the differential Galois methods in a study 
of the integrability. During last ten years we applied them to study several Hamiltonian 



and non-Hamiltonian systems which appear in physics and astronomy, see, e.g., [6; 20:121 
2i S S B S S S 113]. Always a hidden dream of those investigations was a strong 
will to find an unexpectedly integrable system. However, for a long time, neither we, nor 
other authors succeeded with this respect. All those successful investigations gave negative 
results: the investigated systems are not integrable except already known integrable cases. 

In this paper our aim is to give an overview of our works concerning the integrability of 
Hamiltonian systems with homogeneous potentials |24; 2^; 4^ 13; 43], and homogeneous 
Newton equations [41]. Our motivation for those works was an optimistic plan to find the 
necessary and sufficient conditions for the integrability. In other words, our dream was to 
find all possible integrable polynomial potentials and forces. 

Amazingly enough this plan was not only a dream — we found its quite successful reali- 
sation. Our investigations differ in many respects from typical applications of the Morales- 
Ramis theory. In a 'typical' application of this theory there are two difficult problems: the 
first one is to find a particular solution, and the second is to determine the differential 
Galois group of the variational equations along this solution. In this paper we investigate 
'global' multi-parameter problems. We know a priori a class of particular solutions and 
we know also the differential Galois group for a given values of the parameters. The goal 
in our problem is to prove that the system is not integrable for all but finite number of 
parameters' values. To achieve the desired result we developed a method which allows 
to deduce new obstructions for the integrability that come from a global analysis of all 
possible particular solutions of a given class. 

Plan of this paper is the following. In the next section we overview various notions of 
the integrability of ordinary differential equations. In Section [3] we explain how the dif- 
ferential Galois theory can be used for a study of their integrability. The Morales-Ramis 
theory, as well as the Ziglin theory are dedicated for Hamiltonian systems. In our pre- 
sentation we show that, in fact, we can use the differential Galois methods in a general 
context. One can find necessary conditions for the integrability defined adequately to the 
geometry of the considered equations and express these conditions in terms of properties 
of the differential Galois group of variational equations. The next two sections [4] and [5] 
are devoted the integrability analysis of the class of homogeneous potentials. The homo- 
geneity of the potential guarantees the existence of non-equilibrium particular solutions 
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generated by so-called Darboux points. In Section |4] the necessary integrability conditions 
obtained from an application of the Morales-Ramis theory to a particular solution given by 
a Darboux point are formulated. But these conditions are too weak for the ambitious clas- 
sification programme of all integrable potentials. In Section |5] the additional integrability 
conditions obtained from the simultaneous analysis of all Darboux points are formulated. 
In Section |6] it was shown that various parts of analysis made for homogeneous potentials 
can be adapted for systems of Newton equations with homogeneous velocity-independent 
forces. In final Section [7] open problems and perspectives of the classification program of 
integrable systems are discussed. 



2 Integrability 

In this section we overview various notions of the integrability of systems of ordinary 
differential equations. 

Let us consider a system of ordinary differential equations 

^^x = v{x), x={x\...,x^)eR"\ (2.1) 

with smooth right-hand sides v{x) = {v^{x),. . .,v"^{x)). As it was observed a long time 
ago, the knowledge of first integrals or other invariant quantities helps to find explicit 
solutions of this system. Let us recall that a smooth function F is a first integral of system 
jl.l} iff it is constant along its solutions. Thus, a constant value level of F is a set invariant 
with respect to the flow of (|2.1|) . Hence, roughly speaking, knowing a first integral we 
can reduce the dimension of the system by one. Thus, if we know m — 1 first integrals 
Fi, . . . , F„,_i, which are functionally independent in a certain domain U C M"^ then we can 
transform system (|2.1|) into the following one 

^1/ = w{y), w{y) = (0,...,0,it;"'(i/)). (2.2) 

Solutions of this equation can be find easily 

y(t):=yi, for i = l,...,m-\, (2.3) 



and y'"{t) is given by 



y^lO ds 



T = i- (2.4) 

^""{yl yo'^'.s) 

In effect, we reduce the problem to calculation just one integral and inversion of a function. 
In the prescribed situation, we say that the considered system is integrable by quadratures. 
This notion plays a fundamental role. In fact, as we will see later, the other definitions of 
integrability give necessary conditions for the integrability by quadratures. 

A first integral of a given system is just an example of simplest tensor quantity which 
is invariant with respect to the flow of the system. A smooth tensor field F(x) is invariant 
with respect to the flow of the system (|2.Hl iff 

Lv{T) = 0, (2.5) 

where denotes the Lie derivative along vector field v. The existence of a certain number 
of tensor invariants can guarantee that the integration of the system reduces to quadratures. 
Below we give three examples of results of this type. The first of them, due to S. Lie, tells 
that system admitting m linearly independent and commuting symmetries, i.e., invariant 
vector fields, is integrable by quadratures. 
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Theorem 2.1 (S. Lie). Assume that system (|2.1|) admits m independent and commuting symmetries 
Ml = V, U2, ■ ■ ■ , Um- Then the system is integrable by quadratures. 



A differential ra-form co in M'" is given by 




It is invariant with respect to system (|2.1|l iff 



L,(a;)= [f;^g^]dxiA---Adx'" = 0. 




In the classical literature an invariant m-form is called the Jacobi last multiplier. 

Theorem 2.2 (L. Euler, C. G. J. Jacobi). Assume that system (12.11) admits m — 2 functionally 
independent first integrals and invariant differential m-form. Then the system is integrable by 



As we will see later, the following theorem, due to O. I. Bogoyavlensky [10; 11], is very 
important in the context of applications of differential Galois method to the integrability of 
non-Hamiltonian systems. 

Theorem 2.3 (O. I. Bogoyavlensky). Assume that system admits 1 < k < m functionally 
independent first integrals Fi, . . . ,Fj^, and m — k linearly independent and commuting symmetries 
Ml = V, U2, ■ ■ ■ , Um-kf such that 



Then the system is integrable by quadratures. 

Assumptions of the above theorems can serve as a source of definitions for specific 
types of the integrability. Thus, for example, we have the following definition which is 
frequently used in non-holonomic mechanics. 

Definition 2.4. We say that system (|2.H) is integrable in the Jacobi sense iff it admits m — 2 
independent first integrals and invariant differential m-form. 

In a similar way, we can take the assumptions of Theorem l2.3l as a basis for the definition 
of B-integrability. 

Definition 2.5. We say that system {2.1) is B-integrable iff it admits 1 < k < m functionally 
independent first integrals Fi, . . . , Fj., and {m — k) linearly independent and commuting symmetries 
Ml = V,. . . , Um-k, such that My[F,] = 0/or \ <i<k, l<i<m — k. 

The above definition arises from a careful analysis of the concept of the integrability of 
Hamiltonian systems in the Liouville sense (see below). 

To describe shortly the most characteristic features of B integrable systems let us con- 
sider such a system, and let us assume that it admits functionally independent first inte- 
grals Fi, . . . , Fjt- With these integrals we associate the momentum map 



quadratures. 



Proofs and detailed discussion of the above two theorems can be find, e.g., in 



L„ (F,) = M^[F,] = 0, for l<i<k; l<i<m-k. 



(2.7) 




(2.8) 



Let us consider a common level of the first integrals 




(2.9) 
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where / G R*^. If first integrals are independent on Mf, then is a smooth manifold. But, 
even if is connected and compact, its topology can be very compHcated. However, by 
the assumed B-integrabiUty, we know that there exist independent and commuting vector 
fields Ml, ... , Mm_jc tangent to Mf, . Thus, if Mf is connected and compact it is diffeomorphic 
to m — k dimensional torus T'"^^, see Chapter 10, Lemma 2 in fz]. In a neighbourhood of 
Mf we can introduce local coordinates {(p, I), where J G D C M'^, and q) = {fi, ... , (pm-k) 
are angular coordinates on T'"^''^. In these coordinates system {2.1) reads 

= ll = 0. (2.10) 

Thus, a solution of this system has the form 

I{t) := lo cp{t) := iv{Io)t + (p^. (2.11) 

From the above considerations it follows that B-integrability is similar to the integrability 
in the Liouville sense of Hamiltonian systems. In fact, the B-tntegrability was introduced 
as a certain generalisation of the integrability in the Liouville sense. 

Let us assume that system (|2.1|) is Hamiltonian. Then m = In, and v = Xh is a 
Hamiltonian vector field generated by a smooth Hamiltonian function H : M^" M. Here 
we consider M.^" as a linear symplectic space with chosen canonical coordinates x = {q,p), 
where q = {qi, . . ., q„) and p = {pi, . . ., p„). In these coordinates the symplectic form Q is 
following 



n = ^ dqi A dpi. 



i=l 



and the vector field Xh reads 



X -y( — —-— — 
Let us recall the definition of the well known notion. 

Definition 2.6. We say that Hamiltonian vector field Xh is integrable in the Liouville sense iff it 
admits n functionally independent and commuting smooth first integrals Fi, . . . , F„. 

We notice here that a Hamiltonian vector field integrable in the Liouville sense is 
B-integrable. In fact, Xh admits n first integrals, and n = In — n symmetries Xf^, . . . ,Xf„, 
which satisfy 

XF,.[Fy] = {f,vf;} = for 1 < < n, (2.12) 

where { •, •} denotes the Poisson bracket. 

Obviously, a B-integrable Hamiltonian system can be non-integrable in the Liouville 
sense. 



3 General theory 

In this section we show how we can use the differential Galois theory to find necessary 
conditions for the integrability of ordinary differential equations. To deduce such condi- 
tions we have to make several assumptions. The most general one is that the considered 
system as well as the considered first integrals, or other invariants, have 'good' analytical 
properties. Moreover, the theory requires that the 'scalars' form an algebraically closed 
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field. Thus, we assume that this field is just the field of complex numbers C. In effect 
we work with complex functions, complex vector fields, and so on. The above mentioned 
'good' analytical properties mean that the considered tensors are holomorphic at points 
where they are defined. 

Let us consider a complex holomorphic system of ordinary differential equations 

= v(x), xeUC C", t G C, (3.1) 

at 

where U is an open and connected subset of C" . The basic assumption for further consid- 
erations is that we know a particular non-equilibrium solution (p{t) of this system. Usually 
it is not a single-valued function of the complex time t. Thus, we associate with (p{t) a 
Riemann surface T with i as a local coordinate. 

The variational equations along T have the form 

^? = A(0?, A{t) = ^^{cp{t)), ^eTrU. (3.2) 

The entries of matrix A(f) in the above equation are elements of field K := M(r) of func- 
tions meromorphic on T. This field with the differentiation with respect to f as a derivation 
is a differential field. Only constant functions from K have vanishing derivative, so the 
sub-field of constants of X is C. 

It is obvious that solutions of (|3.2|l are not necessarily elements of K"\ The fundamental 
theorem of the differential Galois theory guarantees that there exists a differential field 
L D K such that m linearly independent (over C) solutions of (|3.2|) are contained in L'". 
The smallest differential extension L D K with this property is called the Picard-Vessiot 
extension of K. 

A group 9 of differential automorphisms of L which do not change K is called the 
differential Galois group of equation (|3.2|l . It can be shown that S is a linear algebraic group. 
Thus, in particular, it is a union of a finite number of disjoint connected components. One 
of them, containing the identity, is called the identity component and is denoted by 9°- 

Let ^ = (^^,. . .,^'")^ e L'" he a solution of equation (|3.2|) . and g an element of its 
differential Galois group 9- Then, g{^) := {g{^^),. . .,g(^'"))^ is also its solution. In fact, 
by definition g commutes with the time differentiation, so we have 

^^(?) = g (^^) = g (Mm = g{Mt))g{^) = 

as g does not change elements of K. Thus, if E G GL(n, L) is a fundamental matrix of (|3.2|). 
i.e., its columns are linearly independent solutions of l|3.2|) . then g(E) = EM^, where 
Mg G GL(n, C). In other words, the differential Galois group 9 can be considered as an 
algebraic subgroup of GL(m, C). 

Now, we explain how the existence of first integrals of system 0.1} manifests itself in the 
properties of the differential Galois group of variational equations. At first, we introduce 
some definitions. Let us consider a holomorphic function F defined in a certain connected 
neighbourhood of solution q>{t). In this neighbourhood we have the expansion 

F(<p(0+?) =f/(?)+0(||?||'+i), Fi^O. (3.3) 

Then, the leading term / of f is the lowest order term of the above expansion i.e., /(^) := 
f/(f). Note that /(^) is a homogeneous polynomial of variables f = (^^, . . . , ^'") of degree /. 
If F is a meromorphic function, then it can be written as F = P/Q for certain holomorphic 
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functions P and Q. Then the leading term / of f is defined as / = p/cj, where p and q 
are leading terms of P and Q, respectively. In this case /(^) is a homogeneous rational 
function of ^. 

It is easy to prove that if F is a meromorphic (holomorphic) first integral of equa- 
tion i3.1^ , then its leading term / is a rational (polynomial) first integral of variational 
equation (|3.2|) . If system (|3.1|) has k > 2 functionally independent meromorphic first inte- 
grals fi,. . ./f/c/ then their leading terms can be functionally dependent. However, by the 
Ziglin Lemma we can find k polynomials Gi, . . . ,Gk G C[zi, . . . ,Zk] such that 

leading terms of G;(Fi, . . . , Fj^), for 1 < i < k are functionally independent. 

Additionally, if S C GL(m, C) is the differential Galois group of (|3.2|). and / is its 
rational first integral, then f{g{^)) = /(^) for every g ^ 9, see |j3|:l32l], i.e., / is a rational 
invariant of group S. Thus we have a correspondence between the first integrals of the 
system (|3.1|) and invariants of S- 

Lemma 3.1. If equation 0.1} has k functionally independent first integrals which are meromorphic 
in a connected neighbourhood of a non-equilibrium solution (p{t), then the differential Galois group 
S of the variational equations along cp[t) has k functionally independent rational invariants. 

As it was mentioned above, a differential Galois group is a linear algebraic group, thus, 
in particular, it is a Lie group, and one can consider its Lie algebra. This Lie algebra reflects 
only the properties of the identity component of the group. It is easy to show that if a 
Lie group has an invariant, then also its Lie algebra has an integral. Let us explain what 
the last sentence means. Let g C gl(m, C) denote the Lie algebra of S- Then, an element 
V G C gl(m, C) can be considered as a linear vector field: x ^ V(x) := V ■ x, for x G C". 
We say that / G C(xi,. . .,Xm) is an integral of g, iff Y{f){x) = d/(x) ■ Y{x) = 0, for all 
VG0. 

Proposition 3.2. If fi, . . . , f^ G C(xi,...,x^) are algebraically independent invariants of an 
algebraic group S C GL(m, C), then they are algebraically independent first integrals of the Lie 
algebra g o/S. 

The above facts are the starting points for applications of differential Galois methods to 
a study of the integrability. 

If the considered system is Hamiltonian, then we have additional constrains. First of all, 
the differential Galois group of variational equations is a subgroup of the symplectic group 
Sp(2n, C). Secondly, commutation of first integrals imposed by the Liouville integrability 
implies commutation of first integrals of variational equations. The following lemma plays 
the crucial role and this is why it was called The Key Lemma, see Lemma III.3.7 on page 
72 in 

Lemma 3.3. Assume that Lie algebra q C sp(2n, C) admits n functionally independent and com- 
muting first integrals. Then g is Abelian. 

Hence, if g in the above lemma is the Lie algebra of a Lie group S, then the identity 
component S° of 9 is Abelian. 

Using all these facts Morales and Ramis proved the following theorem ||32l:l33ll . 

Theorem 3.4 (Morales-Ramis). Assume that a Hamiltonian system is meromorphically integrable 
in the Liouville sense in a connected neighbourhood of a phase curve T. Then the identity component 
of the differential Galois group of the variational equations along I is Abelian. 
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If the considered system is not Hamiltonian, then there is no a commonly accepted 
definition of the integrability. However, if we restrict yourself to B-integrability, then we 
have a beautiful generalization of Theorem l3.4[ Namely, with system (|3.ip we consider also 
its cotangent lift, i.e., a Hamiltonian system defined in C^"' by Hamiltonian function 

m 

H = J2yi^\x) (3.4) 

where {x,y) = {x^,. . . , x"\ 1/1, . . . , y,„ ) are canonical coordinates in C^"' . Thus, the Hamilto- 
nian equations have the form 

Let us assume that system (|3.H) is B-integrable with k first integrals Fi, . . . , Fj., and m — k 
symmetries Mi, . . . , u^^-k- Then, we claim that Hamiltonian system (|3.5|) is integrable in the 
Liouville sense. Let us define the following functions 

m 

h+j{x,y) := {y,Uj{x)) := ^i/,Mj(x) for 1 < j < m - k. (3.6) 

/=i 

We have 

{Fj+kM} = E (^a^ - ) = iV' l^iM)- (3.7) 

But, [uj, z;] = for 1 < i < m — k, hy assumption, and so Fy+;c are first integrals. We 
show that first integrals Fi,. . . ,Fm pairwise commute. Obviously, we have {F„ F,} = 0, for 
1 < i, i < k. Moreover, we have 

{Fy+,t/f(+)c} = {y, [uj,Ui]) =0 for 1 <i,i <m- k, (3.8) 

as, by assumption, [uj, Ui] = 0. Finally, we have also 

{Fi, Fj+k} = Uj[Fi] =0, for 1 < i < k, 1 <i<m-k, (3.9) 

because, by assumption, Fj, . . . , Fj- are common first integrals of mi, . . . , Uf„-k- Thus, we 
proved our claim. 

Now, let q>{t) = {(p^{t), . . . , <p'"(0) be a particular solution of (|3.1|) . Then, 

t^^{t) := (<p(f),0) ec^'", 

is a particular solution of the Hamilton equation (|3.5|) . The variational equations for this 
solution have the following form 

= Am = -A(O^V. Mt) = ^(Vit))- (3.10) 

The first of the above equations is just the variational equation (I3.2D , and the second one 
is its adjoint. Thus, if S is a fundamental matrix of the first equation in (|3.10|) , then 
X := (E^^)^ is a fundamental matrix of the second equation in (|3.10|) . In effect, the differ- 
ential Galois group of system (|3.10|) coincides with the differential Galois group of its first 
equation, i.e., with the differential Galois group of the original variational equtions (|3.2|) . 
Using the above facts Ayoul and Zung proved in iQ] the following theorem. 
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Theorem 3.5 (Ayoul-Zhung). Assume that system (|3.1|) is meromorphically B-integrable in a 
connected neighbourhood of a phase curve T. Then the identity component of the differential Galois 
group of the variational equations along F is Abelian. 

Let us underline the importance of this theorem. All results which were obtained on 
the basis of Theorem 13.41 and stating that a given Hamilton system is non-integrable in the 
Liouville sense are, in fact, much stronger - the considered systems are not B-integrable. 

Already in the book [32] a very natural extension of described approach was presented. 
Except the variational equations (|3.2|l along the phase curve T corresponding to the partic- 
ular solution (p{t), we can consider also the higher order variational equations. To derive 
them we consider system (|3.1|l in a neighbourhood of F, where we can write the following 
expansion 



where e is a formal small parameter. Inserting the above expansion into equation ((STJI 
and collecting terms of the same order with respect to e, we obtain the following chain of 
equations 

= ^('^-i)), k = 1,2 (3.11) 

where = 0. For a fixed k, this is a system of A:-th order variational equations, and we 
denote it by VEj^. It is a linear and, for k > 1, non-homogeneous system of equations. 
Nevertheless, there exists an appropriate framework allowing to define the differential 
Galois group 5k of VE/f for any k. Obviously, VEi coincides with (|3.2|) , so Si coincides with 
S- A detailed exposition and proofs the reader will find in [36], where among other things 
the following generalisation of Theorem 13.41 is given. 

Theorem 3.6 (Morales-Ramis-Simo). Assume that a Hamiltonian system is meromorphically 
integrable in the Liouville sense in a connected neighbourhood of a phase curve I. Then the identity 
component of the differential Galois group Sjt ofk-th variational equations VEi^ along F is Abelian, 
for all ken. 

Hence, Theorem 13.41 gives only the first order obstructions for the integr ability. If S° = 
S5' is Abelian, then, having only Theorem I3.4[ we cannot be sure whether the system is 
integrable or not. But knowing the above theorem we can continue our investigations and 
check if the S2 is Abelian. If it is not, the system is not integrable, otherwise we have to 
check if S3 is Abelian. This process we continue up to such k that S° is Abelian for i < k, 
and S^ is not Abelian. If we are able to find such k, then the system is not integrable. 

Here it is worth to mention that it is very hard to determine the differential Galois 
groups Sjc with > 1, or even to decide whether Sj° is Abelian or not. This is why we have 



only a few applications of Theorem l3.6l see |3l; 36]. However, all successful applications of 
this theorem show its real power. 

We can derive higher order variational equations for an arbitrary system (|3.1|l . Thus, 
we can ask if we have a generalisation of Theorem 13.51 similar to that described above for 
Theorem l3.6l In fact, we have such generalisation. 

Theorem 3.7 (Ayoul-Zhung). Assume that system (|3.1|) is meromorphically B-integrable in a 
connected neighbourhood of a phase curve T. Then the identity component S^ of the differential 
Galois group S)c of the k-th variational equations along F is Abelian, for all A: G N. 
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4 Integrability of homogeneous potentials - Morales-Ramis theo- 
rem 

Let us consider Hamiltonian systems with n degrees of freedom given by a natural Hamil- 
tonian function 

ii=\tv} + y{^) (4-1) 

where q = {qi, . . . ,q„) and p = [pi, . . .,pn) are canonical coordinates and momenta, V is 
a homogeneous function of degree /: G := Z \ {0}. We assume just from the beginning 
that the considered system is complex, i.e., the phase space is C^". The Hamilton equations 
have the canonical form 

where V'{q) := grad V{q). Morever, we assume also that the time Hs a complex variable. 

We say that a potential V is integrable iff the Hamiltonian system (|4.2|) is integrable. 

One of the most beautiful applications of the Morales-Ramis Theorem 13.41 concerns 
Hamiltonian systems of the prescribed above form. The basic assumption in this applica- 
tion is that there exists a non-zero vector d G C" such that 

V\d) = d, (4.3) 

It is is called a proper Darboux point of potential V. It defines a two dimensional plane in 
the phase spacs C^", given by 

U{d) := {{q,p) eC^" \ q = (pd, p = xpd, {(p,xp) G } . (4.4) 

This plane is invariant with respect to the system (|4.2|) . Equations (|4.2|) restricted to YI{d) 
have the form of one degree of freedom Hamilton's equations 

^,9 = ^, ^t^ = -f'-'' (4-5) 
with the following phase curves 

■■= { if. e C2 I ^xp^ + 1/ = £ I C C2, £ G C. (4.6) 

In this way, a solution {(p,tp) = {(p{t),ip{t)) of (|4.5|) gives rise a solution {q{t),p{t)) := 
{(pd,ipd) of equations (|4.2|) with the corresponding phase curve 

^k,e := {{q,p) e C^" I iq,p) = {(pd,^d), {cp,xp) G F^,, } C U{d). (4.7) 

Morales and Ramis obtained necessary conditions for the integrability in the Liouville sense 
by an analysis of the variational equations along an arbitrary phase curve 1;^ ^ with e ^ 0. 
These variational equations have the form 

x = -(p{tf-^V"{d)x, (4.8) 

where V"{d) is the Hessian of V calculated at d. Let us assume that V"{d) is diagonalisable. 
Then, without loss of the generality, we can assume that V"{d) is diagonal, and in such a 
case system (|4.8|) splits into a direct product of second order equations 

Xi = -Xi(p{tf-^Xi, \<i<n, (4.9) 

where Ai, . . . , A„ are eigenvalues of V"{d). It is easy to show that d is an eigenvector of 
V" {d) with eigenvalue k — 1. We always denote this eigenvalue as A„. 



In [341] J. J. Morales-Ruiz and J. R Ramis proved the following theorem. 



10 



Theorem 4.1 (Morales-Ramis). Assume that the Hamiltonian system defined by Hamiltonian 
(I4.1I I with a homogeneous potential V S C(^) of degree A: G Z* satisfies the following conditions: 

1. there exists a non-zero d G C" such that V'{d) = d, and 

2. matrix V"{d) is diagonalisable with eigenvalues K\,..., A„; 

3. the system is integrable in the Liouville sense with first integrals which are meromorphic in a 
connected neighbourhood U of phase curve I^^^^ with e 7^ 0, and independent on U\ Fj; £. 

Then each pair {k, for i = 1, . . . ,n belongs to an item of the following list 



case 


k 


A 






1. 


±2 


arbitrary 






2. 


k 


P + ^P(P-1) 






3. 


k 


Z. \ K. J 






4. 


3 




1 

-24 + 










1 

-24 + 




5. 


4 








6. 


5 




9 

-40 + 




7. 


-3 




25 

24 - 










25 

24 - 




8. 


-4 








9. 


-5 




49 

40 - 





(4.10) 



w/zere p /s an integer. 

Let us remark that the above theorem does not give any obstruction for the integrabiUty 
if = 2 or /: = -2. 

Remark 4.2 It was explained in iH/ that the assumption that V"{d) is diagonalisable is irrelevant. 
That is, the necessary conditions for the integrabiUty are the same: if the potential is integrable, then 
each A G spectry"(d) must belong to appropriate items of the above list. Additionally, ifV"{d) 
is not diagonalisable, then new obstacles for the integrabiUty appear. Namely, if the Jordan form of 
V"{d) has a block 

"A 1 0' 



/3(A) := 



A 1 
A 



then the system is not integrable. Moreover, if the Jordan form ofV"{d) has a two dimensional 
block /2(A), and A belongs to the second item of table (|4.10|l , then the system is not integrable. This 
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fact was proved in ll^J. In other words, for k ^ {—2,0, +2}, the presence of a proper Darboux 
point d for which V"{d) has a block of dimension greater than two, or block of dimension two 
with corresponding A, belonging to the second item of table (|4.10|) , implies immediately the non- 
integrability of the potential. 

Remark 4.3 The case of a homogeneous potential of degree k = needs a special treatment. Neces- 
sary conditions for the integrability in this case were found recently in ill/ . 

We denote by Mjt a subset of rational numbers A specified by the table in the above 
theorem for a given k, e.g., for \k\ > 5, we have 



Mic = + I pezjuj^ ^ + p{p + i)k I pezj 



(4.11) 



The Morales-Ramis Theorem 14. II is a powerful tool. In fact, the necessary conditions for 
the integrability of polynomial potentials are reduced to solving algebraic equations: we 
have to find a Darboux point and then to check if the eigenvalues of the Hessian at this 
Darboux point are rational numbers which belong to the Morales-Ramis table. However, 
we note that none of these algebraic problems is trivial. First of all, even if a Darboux point 
exists, generally we are not able to find its coordinates. Moreover, we have much more 
serious problems in a parametric case which is the most important in applications. This is 
illustrated by the following example. 

Let us consider the following potential 

V = l^ql + ^qlq2 + ^cql (4.12) 

where a and c are in general complex parameters. For generic values of these parameters 
V has three Darboux points di, di, and dj. The non-trivial eigenvalues A, = TrV" [di) — 2 
of Hessian V"{q) at these Darboux points are following 

1 2c - 1 2c - 1 

Ai = -, A2 = 7— 3— r, A3 = — -, (4.13) 

c 1 + fl^ + A 1 + fl^ — A 

where 



A = -^fl2(2 + fl2_2c). 
If potential V is integrable, then by Theorem 14. 11 we have 

6 



Ai e M3 = U M^'\ for 1 < z < 3, (4.14) 



where 



mW := |p + ^p(p-l) I pezj, M(2):=|10+3p(p + l)) | pGzj, 
{^(l + 3p)'-^| pez}, M(4):=||(l + 4p)2-l| pezj 



From Eq. I|4.13|) we find that 

1 Ai + AiA,- 



Ai V2AiA,(2-Ai-A,) 



for /e{2,3}. 
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Hence, for arbitrary A-[, A2, A3 G M3, the above defined values of a and c give potential (|4.12|) 
which satisfies the necessary conditions for the integrability. 

Theorem 14.11 gives only necessary conditions for the integrability, and there are many 
examples that they are not sufficient. Thus in the above example, we have to check whether 
infinitely many potentials are integrable or not. 

5 Integrability of homogeneous polynomial potentials. Global 
analysis 

In this section we assume that the considered potential V{q) is polynomial and homoge- 
neous of degree k > 2. The set of homogeneous polynomials in n variables q = {qi, . . . ,q„) 
of degree k, we denote by Cjt [q] . 

During last few years we worked on the following problem. Is it possible, for a given 
k > 2 and n > 2, to distinguish all meromorphically integrable potentials V G C|f[^]? In 
other words, is it possible to give a necessary and sufficient conditions for the integrability 
of homogeneous polynomial potentials? The example given in the end of the previous 
section shows that, except Theorem 14. 1[ we need a result which gives stronger necessary 
conditions. This example shows that even for fixed k and n. Theorem 14.11 distinguishes 
infinitely many parameters' values for which the potential can be integrable. 

5.1 Darboux points 

It is clear that the more Darboux points of given potential we know, the more obstructions 
for its integrability we obtain from the Morales-Ramis Theorem 14.11 Hence, we have to 
know how many Darboux points a polynomial potential of a given degree can have. An 
analysis of this and similar problems related to particular solutions of Hamiltonian systems 
with homogeneous potentials forced us to give a more geometrical definition of Darboux 
points. 

Let y be a homogeneous polynomial potential of degree k > 2, i.e., V ^ Cj^lq]. A 
direction, i.e., a non-zero d G C", is called a Darboux point of V iff the gradient V'{d) of V 
at d is parallel to d. Hence, d is a Darboux point of V iff 

dAV'{d)=0, d^O, (5.1) 

or 

V'{d) = jd, d^ 0, (5.2) 

for a certain 7 G C. Obviously, if d satisfies one of the above conditions, then d = ad 
for any a G C* satisfies them. However, we do not want to distinguish between d and d. 
Hence we consider a Darboux point d = {di, . . . , d„) G C" as a point [d] := [di : ■ ■ ■ : d„] in 
the projective space CP"^^. 

The set Ti{V) C CP"^^ of all Darboux points of a potential y is a projective algebraic 
set. In fact, D{V) is the zero locus in CP"^^ of homogeneous polynomials G C^iq] 
which are components of </ A V'{q), i.e. 

dV dV 

Rii := qi- t?/^, where 1 < z < / < n. (5.3) 

oqj ' dqi 

We say that a Darboux point [d] G Ti{V) is a proper Darboux point of V, iff V'{d) ^ 0. The 
set of all proper Darboux points of V is denoted by D*{V). If [d] eD{V)\D*{V), then [d] 
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is called an improper Darboux point of potential V. We say that [d] is an isotropic Darboux 
point, iff 

4 + • ■ ■ + = 0. (5.4) 

We say that potential V is generic iff all its Darboux points are proper and simple. The 
basic fact concerning Darboux points of generic potentials is given in the following lemma. 

Lemma 5.1. The set of generic potentials Qn,k C C^lq] of degree k is a non-empty open set in C;t[^]- 
A generic V [q] has 

proper Darboux points. 

A non-generic potential can have finite, or infinite number of Darboux points, but for 
an arbitrary V G Ck[q] the set Ti{V) is not empty. Moreover, if V does not have improper 
Darboux points, then it has a finite number of proper Darboux points. 



5.2 Obstruction for the integrability due to improper Darboux point 

Here we must justify the introduced definition of the Darboux point. Let us notice that 
in Theorem 14.11 only proper Darboux points appear and they give particular solutions. 
However, we have a more general fact. 

Lemma 5.2. If [d] is a proper Darboux point of a homogeneous potential V of degree k > 2, then 

q{t) := (p{t)d, p{t) := (p{t)d, (5.5) 

is a solution of Hamilton's equation (|4.2|l provided ip = —cp^^^. Moreover, V"{d) ■ d = A„rf with 
An = k — 1, and if additionally [d] is isotropic, then A„ is a multiple eigenvalue ofV"{d). 
If [d] is an improper Darboux point, then 

q{t) := td, p{t) := d, (5.6) 

is a solution of Hamilton's equations (|4.2|) . Moreover, V"{d) ■ d = \nd, with A„ = 0, and if 
additionally [d] is isotropic, then A„ is a multiple eigenvalue ofV"{d). 

Hence, also an improper Darboux point gives the particular solution i5.6^ of the consid- 
ered canonical equations (|4.2|) . However, this solution has an extremely simple form and 
one can doubt if using it we can obtain any obstruction for the integrability. In fact, it is 
easy to notice that the monodromy group of the variational equations along solution (|5.6|l 
is trivial. Thus, in the frame of the Ziglin theory we do not obtain any obstacle for the 
integrability. Nevertheless, in l|42ll the following theorem was proved. 

Theorem 5.3. Assume that a homogeneous potential V G C^lq] of degree k > 2 admits an improper 
Darboux point [d] G CP"^^. IfV is integrable with rational first integrals, then matrix V"(d) is 
nilpotent, i.e., all its eigenvalues vanish. 

For n = 2, the above theorem coincides with Theorem 2.4 in [27]. 
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5.3 Relation among eigenvalues 

For a Darboux point [d] G "^{V) we can calculate eigenvalues Ai(rf), . . . , A„(rf) of the 
Hessian matrix V"{d). However, numbers A,(d) are not well defined functions of point 
[d] e CP"^^, as they depend on its representative d. There are several possibilities to de- 
fine properly the quantities related to the eigenvalues of V"{d) which do not depend on 
a choice of a representative of [d]. However, because of some historical reasons and the 
convention widely accepted in the literature, we choose the one which is a simple normal- 
isation. Namely, if [d] is a proper Darboux point, then the chosen representative is such 
that it satisfies V'{d) = d. If [d] is an improper Darboux point, then the representative of 
[d] can be chosen arbitrarily. 

Let [d] be a proper Darboux point of potential V. Then, thanks to our assumption, the 
eigenvalues Ai{d),. . .,\n{d) of the Hessian matrix V"{d) can be considered as functions 
of d. According to our convention A„(fi) = /c — 1 is the trivial eigenvalue. Let X{d) = 
{Xi{d),. . . , A„_i [d) ) . Hence we have the following mapping 

3 [d] I — > X{d) e C"-\ (5.7) 

Assume that 'D*{V) is finite. Then the image of T)*{V) under the above map is a finite 
subset of C"^^. The question is if we can find a potential V of degree k such that the 
elements in the image have values prescribed in advance. We show that the answer to 
this question is negative. More precisely, we prove that among \{d) taken at all proper 
Darboux points [d] G T)*{V) a certain number of universal relations exists. These relations 
play the central role in our considerations. 

To formulate our first theorem we define A{d) = {Ai{d), . . . ,A„_i{d)), where A,(d) := 
\i{d) — 1 for / = 1, . . . , n — 1. By T,. for < r < n — 1, we denote the elementary symmetric 
polynomials in (n — 1) variables of degree r, i.e., 

r 

Tr{x) ■.= Tr{Xi,...,Xn-l) = ^ Yl^'^' 1 < T < Tl - 1, 

l<ii<---<ir<n-ls=l 

and To(x) := 1. 

Our first theorem gives the explicit form of the above mentioned relations among A{d), 
for a generic potential V. 

Theorem 5.4. Let V ^ Ck [q] be a homogeneous potential of degree k > 2 and let all its Darboux 
points be proper and simple. Then 

and 

for r = Q, . . . ,n — 1. 

Let us explain the importance of Theorem |5.4[ To do this we need more definitions. 

Let C„, denote the set of all unordered tuples A = (Aj, . . . , Am), where A, G C for i = 
1, . . . ,m. For M > 0, the symbol denotes the set of all unordered tuples ( Ai, . . . , A^), 
where A, G Qm, for z = 1, . . . , M. 
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We fix A; > 2 and n > 2, and say that a tuple A G Cn-i is admissible iff A, = A,- + 1 G M)t 
for z = 1, . . . , n — 1. In other words. A, is admissible iff A,- + 1 belongs to items, appropriate 
for a given k, in the table of the Morales-Ramis Theorem 14. 1[ for / = l,...,n — 1. We 
denote the set of all admissible tuples by A„^]^. If the potential V is integrable, then for 
each [d] G T)*{V), the tuple A{d) is admissible. The set of all admissible elements Anjc is 
countable but infinite. 

If the set of proper Darboux points of a potential V is non-empty, and N = card D*(y), 
then the N-tuples 

Ji{V) := {A{d) I [d] G D*{V) ) G e~ i, (5.10) 

is called the spectrum ofV. Let A^^. be the subset of consisting of N-tuples (Ai, . . . , Ajv), 
such that A, is admissible, i.e.. A, G ^n,*:/ for / = 1, . . . ,N. We say that the spectrum ^{V) 
of a potential V is admissible iff ^{V) G -AJ^^^. The Morales-Ramis Theorem 14.11 says that 
if potential V is integrable, then its spectrum ( F) is admissible. However, the problem is 
that the set of admissible spectra A^^. is infinite. We show that from Theorem I5.4l it follows 
that, in fact, if V is integrable, then its spectrum ^{V) belongs to a certain finite subset Jj^j, 
of A^j^. We call this set distinguished one, and its elements distinguished spectra. 

Theorem 5.5. Let potential V satisfy assumptions of Theorem 15.41 If V is integrable, then there 
exists a finite subset jj^^ C A^^, where N = card'D*{V), such that Si{V) G 3^j^. 

Informally speaking, for fixed k and n, we restrict the infinite number of possibilities in 
each line of the Morales-Ramis table to a finite set of choices. 



5.4 Euler-Jacobi-Kroncker formula and its generalisation 

The importance of Theorem 15.51 is clear. Having in mind the our general program of find- 
ing all integrable potentials, one would like to have a generalisations of Theorem 15.41 and 
Theorem 15.51 for non-generic potentials. In [40] we gave a proof of Theorem 15.41 using a 



certain result of Guillot [17]. Unfortunately, the methods used in |17f ] do not admit such a 
generalisation. 



Our analysis of case n = 2 given in [27l ] explicitly showed that one can find an alterna- 
tive proof of Theorem 15.41 which admits a generalisation to non-generic cases. Moreover it 
also gives a clue that an alternative proof Theorem 15.41 can be done with a help of multidi- 
mensional residue technique. We have made many attempts to find such a proof, however 
all of them failed. 

Finally we have found amazingly simple solution of the problem. Here we describe 
shortly the main construction of our approach. A detailed exposition the reader will find 
in Blii. 



Let us introduce local affine coordinates on CP" ^ where the Darboux points live. We 
choose chart (Ui, where 

Ui := CP"-i \ { [^] e Cr'-i I qi^o}, 

and 

0i:!Ji^C"-i, x:={x,,...,x„^,) = 9,{[q]), (5.11) 

where 

Xi = ^, for / = l,...,n-l. (5.12) 
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The image of the set of Darboux points which lie on this chart, i.e., 0i(D(y) fl Ui), is an 
affine algebraic set 

9i{'D{V)nU^)=V{gi,...,gn-i), (5.13) 
where polynomials gi, . . ■ ,gn-i S C[x] are given by 

v{x) ■.= ¥{!, xi,...,Xn-i), gQ := kv - Yj > (5-14) 

i=i 

and 

dv 

gi •= 9^ - ^'^0, for z = 1, . . . , n - 1. (5.15) 

Moreover, [d] G D(y) fl Ui is an improper Darboux point iff its affine coordinates a := 
Biild]) satisfy ^o(«) =0. 

It is instructive to consider first the case n = 2. In this case a Darboux point [d] on the 
affine chart is given by one coordinate = 0i ( [d] ) . It is a root of polynomial 

gi{x) := v'{x) — xgo{x) where goi^) '■= kv{x) — xv'{x). 

Moreover, it is easy to notice that the non-trivial eigenvalue A((i) of V"{d) can be calculated 
from the following formula 

1 - ^ -S'^'*^ (5.16) 



A{d)-1 A{d) g[{x 
The above formula suggests to introduce the following differential form 

a; = 44dx, 

considered as a differential form on CP^. This form has poles at Darboux points. If a 
Darboux point [d] is proper and simple, then its affine coordinate x+ is a simple pole of co, 
and the residue of co at this point is 

go{x^) 1 
res(a;, x+) — ° ' ' — 



Without loss of the generality we can assume that all Darboux points are located in the 
affine part of CP^, and then 

res(a;, oo) = 1. 

Thus, assuming that all Darboux points are proper and simple and applying the global 
residue theorem we obtain that 

E ^ = 1- 



This is just relation (|5.8|) for n = 2 with r = 0. Note that for n = 2 it is the only non-trivial 
relation. 

Now, considering cases with n > 2 one would like to construct differential forms in 
CP"^^ which have poles at Darboux points and such that their multidimensional residues at 
these poles are given by symmetric functions of A(d). It is not difficult to define an (n — 1)- 
differential form in affine part of CP"^^ which has poles at Darboux points with residues 
given in term of symmetric functions of A{d). The problem appears with extension of 
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this form onto whole CP"^^. The obtained global differential form has some additional, 
usually not isolated poles and for this type of forms there is no an appropriate global 
residue theorem. 

Let us recall basic facts about the multi-dimensional residues and the Euler-Jacobi- 
Kronecker formula. For details the reader is refered to Q; [H; E H Si. Let /, : C" D 
JJ — > C, where U is an open neighbourhood of the origin, be holomorphic functions for 
z = 1, . . . , n, and x = be an isolated common zero of We consider differential n-form 

p(x) 

CO := , , , ' , , dxi A ■ ■ ■ A dx,,, (5.17) 
/i(x) • ■■jn{x) 

where p : iJ C is a holomorphic function. The residue of the form a; at x = can be 
defined as 

where 

r:={xG!J| |/i(x)| =£i,...,|/„(x)| =£„}, (5.19) 
and £i, . . . ,£„ are sufficiently small positive numbers. The orientation of F is fixed by 

d(arg/i) A ■ ■ ■ A d(arg/,) > 0. (5.20) 

Let us denote /:= (/i, ...,/„). It can be shown that if the Jacobian det/'(0) 7^ 0, then 

The following theorem gives the classical Euler-Jacobi-Kronecker formula, see e.g. H^. 

Theorem 5.6 (Euler-Jacobi-Kronecker). Let fi,...,fn £ C[x] be non-constant polynomials such 
that V(/) := V{fi, . . . ,fn) is finite and all points of this set are simple. Iffi, . . . ,fn do not intersect 
at the infinity, then for each p S C[x] such that 

degp<£deg/,-(n + l), (5.22) 

(=1 

we have 

dGV{f) dGVif) ^^^J V"J 

The above theorem is not sufficient for our investigations. We have to consider cases 
when fi,...,fn have intersections at the infinity as well as cases when intersections of 
fi,...,f„ are not simple. 

The homogenisations of // are given by 

F,(zo,Zi,. . .,z„) := Zq"^-^'/,- (—,...,— ) , for i = l,...,n. (5.24) 

They define the projective algebraic set V{F) := V{Fi, . . . ,F„) C CP" whose affine part 
is homeomorphic to Next we extend the form a; to a rational form Q defined on 

CP". To this end we consider co as the expression of O on the chart {Uo,6o). In order to 
express O on other charts we use the standard coordinate transformations of n-form. Let 
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[p] = [po '■ ■ ■ ■ '■ Pn] £ Ui n V(fi, . . . , F„). We can define the residue of the form O at point 

[p] as 

res(a[p]) :=res(c5,0,(M)), (5.25) 

where co denotes form Q expressed in the chart (U,, 0,). 

The form Q is defined by homogeneous polynomials fi, . . . , f,„ and 

P(zo,Zi z„):=zJ^SPp/'£i (526) 

To underline the explicit dependence of Q on f, and P we write symbolically O = P/F. 
The following theorem is a special version of the global residue theorem. 

Theorem 5.7. Let V{F) := V(fi, . . . , F„) be a finite set. Then for each polynomial P such that 

degP< £degF, -(n + l), (5.27) 

i=l 

we have 

X] res(P/f ,[s]) =0. (5.28) 

Wgv(f) 

For the proof and the more detailed exposition we refer the reader to |@; 15]. 

If G "^{f) is an isolated but not simple point, then we cannot use formula {5.11} to 
calculate the residue of the form to at this point. In such a case we can apply a very nice 
method developed by Biernat in i't'; that reduces the calculation of multi-dimensional 
residue to a one dimensional case. We describe it shortly below. 

Let us consider the following analytic set 

A:={xeU \ f2{x) = ■■■= fn{x) =0}, (5.29) 

where U C C" is a neighbourhood of the origin. Set A is a sum of irreducible one di- 
mensional components ^1 = ^li U • • • U Am- Let t ^ q)-{t) S Ai, ^,(0) = 0, be an injective 
parametrisation of Ai. Then we define the following forms 

/(?>,■(<)) fM(t)) 

As it was shown in [8] we have 

m 

res(a;,0) = ^res(a;,,0). (5.31) 

i=l 

In order to use the above theorems we have to make a kind of blowup. Roughly speak- 
ing, the idea is to associate with a Darboux point which is located in CP"^^, a finite set of 
points in CP". 

Let us define the following n homogeneous polynomials of n + 1 variables q := {qo, qi, . . . , 

f; := — - q^-^qi, i = l,...,n, (5.32) 
dqi 



and an algebraic set T>{V) = V(fi, . . . , f„) C 

Assume that [d] G D*(y). Then there exists 7 G C*, such that V'{d) = jd, so k - 2 
points [ '^^^ : di : • • ■ : d„] S CP" belong to D(y). These points are well defined as they do 
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not depend on a representative for [d] . If [d] is an improper Darboux point, then it defines 
just one point [0 : di : ■ ■ ■ : d„] S CP" which is a point of T>{V). 

Set is not empty because it contains point [do] := [1 : : ■ • • : 0]. If [d] = [do : ^^i : 

• • • : d„] S 23(y) \ {[do]}/ then [d] = [di : ■ ■ ■ : d„] is a Darboux point of V. Moreover, if 
do 7^ 0, then [d] is a proper Darboux point. 

The natural projection 



7r:CP"\{[do]} 



n-l 



TT{[qo : (ji : ■ ■ ■ : qn]) = [^i ■ ■ ■ ■ ■ In], 



(5.33) 



maps T){V) \ {[do]} onto D{V), and the intersection of the inverse image n ^{[d]) of a 
Darboux point [d] G with T>{V) is a finite set. We define also 



n:%{V)\{[do]}^'D{V), 



(5.34) 



putting 7r([(i]) := /T([d]) for [d] G I'(y) \ {[t^o]}- That is, n is the restriction of tz to 
'D{V) \ {[do]}- This construction is illustrated in the Figure [H Now, we can consider 





Figure 1: Sets C CP" and V{V) C CF"-'^ 

differential form Cl := P/F in CP", which in affine part of CP", is given by 

Pr{x) 



^' fl{x) ■ ■ ■ fn{x)' 

where /, is a dehomogenisation of F,, i.e., 

dV 

fi{x) := — (x) - Xi, for i = l,...,n, (5.35) 

and polynomials pr are of the form 

Pr{x) = {Tif{x) - {k-2)y, with r e {0,...,n-l}. (5.36) 

Now, to obtain relations (|5.8|) we just make simple calculations according to the following 
scheme 

1. V(/) :=V(/i,...,/„) is the affine part of V(F) := V(Fi, . . . ,F„) = S(y). 
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2. G V{f), and 



f{0) = -En so: pr{0) = {-n-{k-2)y, det/(0) = (-1)". 

3. If d G V{f), and d 0, then [d] G V*{V) and 

det/(d) = (fc-2)nA,(d), and = ^"^'A,(d)^ . 

4. For k > 2, and n >2 we have 

n 

detpr = r{k-2) < X]deg/, - (n + 1) = n(/: - 1) - n - 1, 
for r G {0,. . .,n - 1}. 

5. Moreover, := t^d G V(/), and f {dj) = f{d) for ; = 0, . . . ,k — 3, where e is a 
primitive {k — 2)-root of unity. 

6. If all Darboux points are proper, then polynomials fi,. . . ,fn do not intersect at the 
infinity and we can apply the classical Euler-Jacobi-Kronecker formula (|5.23D from 
Theorem |5^ 

In a case of non-generic potential V having finite number of Darboux points the potential 
possesses either improper, or multiple Darboux points. In such cases we apply Theorem l5.7[ 
Examples of such calculations are given in [43t]. 

5.5 Applications of global analysis 

In order to perform a reasonable classification of potentials it is convenient to introduce the 
following equivalent relations. 

Let PO(n, C) be the complex projective orthogonal subgroup of GL(n, C), i.e., 

PO(n,C) = {A G GL(n,C), | AA^ = ocE„, a G C*}, (5.37) 

where is n-dimensional identity matrix. We say that V and V are equivalent if there 
exists A G PO(n,C) such that V{q) = Va(<7) := V{Aq). Later a potential means a class of 
equivalent potentials in the above sense. 

The general results described in the previous section can be applied to a systematic 
study the integrability of homogeneous potentials with fixed n and k. The algorithm is 
following. 

We assume that n > 2 and k > 2 are fixed. The aim is to distinguish all integrable 
potentials. 

At first we consider generic potentials with N = D{n,k) proper Darboux points. 

By Theorem 15.51 there is only a finite number of distinguished spectra J^j, and we 
can find all them solving Diophantine equations of the form (|5.8|l in a subset of rational 
numbers defined by the Morales-Ramis table. For example, for n = 2 we have D{2,k) = k, 
so a generic homogeneous potential of degree k, has k proper Darboux points. At each 
proper Darboux point we have one non- trivial eigenvalue A. Thus, in this case elements 
of a Jj^j, = j\ ^, are unordered tuples of k elements. For k = 3 and k = 4 they are listed in 
Table [1] 
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(-1,-1,1) 
(-2/3,4,4) 
(-7/8, 14,14) 
(-2/3,7/3,14) 



(Ai,A2,A3,A4) 

(-1,-1,2,2) 
(-5/8,5,5,5) 
(-5/8,2,20,20) 
{-5/S,27/S,27/S, 135) 
(-5/8,2,14,35) 



Table 1: The distinguished spectra "J^^ for k = 3 (left), and A: = 4 (right) 

It appears that the determination of the distinguished spectra can be performed only 
with a help of a computer algebra system. The known algorithms used for these purposes 
are highly time demanding. 

The next step is to find all possible potentials for a given element of the distinguished 
spectrum. In other words we have to determine all polynomials of a fixed degree such that 
their Hessians at some points (which are unknown a priori) have specified eigenvalues. At 
a first glance, it seems that this problem is ill-posed. However, it is not like that. The reason 
is that the restriction imposed by the fixing of the eigenvalues is very rigid. Moreover, in 
fact we work only with equivalent classes of potentials and this restricts additionally the 
number of free parameters in the problem. The algorithm of performing this step is based 
on determination of the elimination ideal. In all cases, for a fixed distinguished spectrum 
we obtained, either a finite number of non-equivalent potentials, or a finite number of 
families of 'separated' potentials (a sum of potentials which depend on smaller number of 
variables). 

From the previous step we obtained a finite number of potentials which can be inte- 
grable. To check if they are integrable we use two methods. First we try to find an addi- 
tional polynomial first integral applying a direct method. If it fails we apply the higher 
order variational equations in order to prove their non-integrability. 

In a similar way we investigate non-generic cases. 

The first time the prescribed algorithm was used in | |23 ] where it was shown that all 
integrable homogeneous of degree three polynomial potentials in two variables are already 
known. Next, in [27] it was shown for n = 2 and /c = 4 all integrable potentials are known 
except for potential 

with OL such that A = 1 — a belongs to items 2,3 and 5 of table (|4.10l >. All those investigations 
of homogeneous potentials with two degrees of freedom allowed to find one new non- 
trivially integrable potential of the form 

Vk,i = {q2-m)\q2 + iqif'K (5.38) 

with k = 7 and I = 2, which admits an additional polynomial first integral of degree four 
in the momenta, see [38]. 

The case of homogeneous degree three polynomial potentials in three variables, i.e.. 



case n = k = 3, was analysed in ||42l: |43[|. In this case a generic potential has seven proper 
Darboux points [d;]. At each of them we have pair A{d,) := (Af'^A^')) of the shifted 
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eigenvalues. The relations (|5.8|) have the following form 



E 



1 aS')a« 



1, 



E 



7 aS') + a 



(0 

2 



- aS')aJ-) 

f (Af)+A^ 



-4, 



16. 



(5.39) 



After really long computations we have found ten distinguished spectra. For each of them 
we reconstructed the general form of the potential. Four among the reconstructed poten- 
tials have the form 

1 , 

where Vi{qi, q2) is an integrable potential with two degrees of freedom. Each distinguished 
spectrum gives several potentials. Some of these potentials are not integrable, and this fact 
was proved with the help of higher order variational equations. The remaining ones have 
the forms 



iql + 6552iVl7q] 



Vb =^qiqi + -^qi + i^qiq^ + ^qi 

Ve =364:VV7ql + 2835iVV7qlq2 + 1560VV7qic 
4335qlq3 + 19074ij|(?3 + 578ql, 

- + 2mVUqlq2 + 330V7qiql + 935iVUql + 3087qlq3 + 294^/^, 
51^3 2 9iV3 . 15 

— —4' 



Vy =uV7ql 

7 , 
^8 =2^i'?3 



^8 =2^i'?3 — ^qiqi — ^qi + y'?2'?3 + ^qi 

V9 =27W3990ql + 3726v/T5(??(?2 - ^56iV3990qiql - m2VT5q] 
lUSqlqs - 3000qlq3 - 50ql, 

V.0 ' 



3000qlq3 - 50ql, 



All the above potentials are integrable. Each of them admits two commuting additional 
polynomial first integrals Ji and J2. They were found with the help of a direct method. 
All these potentials are integrable in a non-trivial way, i.e., at least one of additional first 
integral is of degree higher than two with respect the momenta. For example, for the 
potential Vio the additional first integrals have the forms 



h = 12pf - 27q^2 - IS'jfl'?? - 4v^(?i(?3 + 2ql) 
+ 12ql{3pl{V2q^ - ^3) + Upipsiqi + V2q3 
-12p2'?2(2p3(16t?? 
-12pi(2p3(2v^pi 



-4(6p? 
-2ql{12ql 



5PI 



+ 16V2ql-2ql){3pl 
V2q^q3 + 2ql)) 



2ql) 



3ql 



8^2qlq3-^l)+3^/2pl{ql+^l)) 



Ps) -4((?2-'?3)'?3('?2 + '?3) - V2qi{Sql + Sql)), 
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h = SlqlilVlq^ + qs) + 216p2P3^?2(v^'?i + 2(?3) + 54rjf - 3p| + 4^/2^?? - lAqlqs 

- eVlqiql) + 384p2P3'??'?2(3pi + sVlql + 8i??i?3 - ^Vlq^qj) - 72p\{3pl + 2ql) 

+ lUp2p3ql{pl + Sq\{lVlqi + 3^?3)) + UApl{V2plp3 + 3V2p2^j2<?3 " 3p3'?i('?i + v^'?3)) 

- 32(pf + 12p4^?<?3 + UplqliVlpl + 4^?i<?i) + 32<??(3pi + 2ql)) - 12pf (4pf 

- 6p2P3q2{16ql + 9^?^ + 8\/2iji(?3 - 4ql) + 9(jf (2ij? + ^Vlqiq^ + + 3lV2q\{3pl + 2(?^) 
+ 12p|(pi + 4.qlq3 - V2qi{ql - 2ql)) + 6ql{9V2plqi + 2^ji(-6^?? + iVlqiqs + ql))) 

- 144(?4(p2(7^2 ^ 5^^^^3 + 2rj|) + 3(?f (3pi - 2q3{-2ql + 2V2^?i(?3 + ^^s))) 

- 48ql{pj{5V2qi + 4^?3) + 4.plql{8ql + 2v^^?i(?3 + 3(ji) + 8ql{9plqi + (?i(-6v^(j? 

+ 4(^i(?3 + Vlql))) + 6piil6V2p^p3 + 16p|p3(8(?? -6(ji(?| + 3v^(j^^?3) +4pi^?2(-16v/2<?5 

+ 32qiq3 + V2{3ql + 4^j§)) + 3p3ql{9^2q\ - 6Aql{V2qi + 2q3) - 12ql{2\/2ql + 8(?iij3 

+ \/2ql)) + 12p2q2{-3V2plql + 9qiq\ + 32rj^i?§ + 4i?^(4i?f + 6q^ql + \/2(?^))). 

For the potentials Vj and V9 integrals 7i and 72 have, as above, degrees four and six with 
respect to the momenta, respectively. However, they are 'monster' first integral — few pages 
are necessary to write them down, see (i?]. 

The analysis of non-generic cases is much more involved but, nevertheless it can be 
made almost up to the end, see |,43. ] for details. 

6 Integrability of Newton homogeneous equations 

In this section we consider the following class of Newton's equations 

q = -F{q), q = [q]^^...^q„)^ (6.1) 

that we rewrite as a system of first order differential equations 

q = p, p = -F{q). (6.2) 

Our aim is to investigate integrability properties of such equations. Generally, the forces 
F are not potential, so it can happen that these equations do not admit even a single first 
integral. 

It is easy to observe that equations (|6.2|) admit 

ji = dqi A • • • A dqn A dpi A • • • A dp,,, 

as an invariant 2n-form. Thus, we can talk about the integrability in the Jacobi sense of 
such equations, see Definition 12.41 

Now, we want to find necessary conditions for the integrability in the Jacobi sense 
applying the differential Galois framework. Thus, we assume just from the beginning that 
{(J,p) S C^", and we consider only the case when forces are homogeneous of degree (k—l). 

Let us notice formal similarities between Newton's (|6.2|l , and Hamilton's (|4.2|l equations. 
These similarities allow to define the notion of Darboux point of a homogeneous force F, 
as a non-zero direction d such that F{d) is parallel to d. We say that <i is a proper Darboux 
point iff F{d) 7^ 0. As in the case of homogeneous potentials, a Darboux point d defines an 
invariant two dimensional plane 

n{d) := {{q,p) eC^" \ q = (pd, p = xpd, {(p,xp) G } . (6.3) 
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Equations (|6.2|) restricted to Yl{d) have the form of one degree of freedom Hamilton's 
equations 

q> = tp, ^J = -cp^-\ (6.4) 



dr ^' dt 
with the following phase curves 



(6.5) 



In this way, a solution {cp,tp) = {(p{t),ip{t)) of l|6.4|) gives rise a solution {(j{t),p{t)) := 
(^rf, t/;d) of equations ( I6.2D with the corresponding phase curve 

:= K'/.P) e C^" I = ((pd,i/;d), {cp,xp) G r^,, } C U{d). (6.6) 

The variational equations along r;t,£ have the form 

x = y, y = -(p{tf-^P{d)x, (6.7) 

or simply 

x = -(p{tf-^P{d)x, (6.8) 

where f '(d) is the Jacobi matrix of F calculated at a Darboux point d. Let us assume that 
this matrix is diagonalisable. Then, in an appropriate basis equations (|6.8|) have the form 

iii = -\i(p{tf-\, i = l,...,n, (6.9) 

where Ai, . . . , A„ are eigenvalues of F'{d). It is easy to show using the Euler identity that 
d is an eigenvector of F' {d) with eigenvalue k — \. We always denote it by A„. In [41] the 
following theorem was proved. 

Theorem 6.1. Assume that the Newton system \6.1) with homogeneous right-hand sides of degree 
{k—l) with k G Z*, satisfies the following conditions: 

1. force F admits a proper Darboux point d, and Ai, . . . , A„ are eigenvalues ofF'{d), 

2. equations (|6.2|) are integrable in the Jacobi sense with first integrals which are meromorphic 
in a connected neighbourhood U of phase curve I^^i^ with £ 7^ 0, and independent onU \ Tk^^- 

Then pairs {k, Aj) for i = 1, . . . ,n belong to the Morales-Ramis Table (|4.10|) . 

The above theorem is a generalisation of two theorems (Theorem 1.2 and Theorem 1.3) 
from fZl*]. Here we remark that the similarity of the theses of Theorem l4.1l and Theorem l6.1l 
is somewhat misleading. The point is that in the case of Newton equations we do not have 
any symplectic structure, and thus we have no all geometrical consequences of this fact. On 
the other hand, the integrability in the Jacobi sense requires 'big' number of first integrals 
and exactly this requirement is the reason why the identity component of the differential 
Galois group must be Abelian. 

Fact that we have Theorem 16.11 allows us to think about a global analysis similar to that 
we presented in Section|5]for the homogeneous potentials. In order to analyse this question 
we assume the force F has polynomial components of the same degree Z := — 1 > 1, i.e., 
F G (C/ [</])". We denote by D(F) the set of all Darboux points of f. It appears that a generic 
force has D{n,k) = [{k — 1)" — l]/{k — 2) Darboux points (considered as points in CP"~^), 
i.e., as many as a homogeneous potential V of degree k. It is rather amazing, because 
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Fv{q) ■= V'{q) depends on ("^^7^) parameters while F depends on n("^^2^) parameters, 
so in the second case the number of parameters is much bigger. 

For [d] S we denote eigenvalues of F'{d) by M{d), . . .,A„{d) = {k — 1), and we 

set A{d) := (Ai(d), . . . , A„_i((i)), where A,(d) := A,(rf) - 1, for i = l,...,n-l. 

For the considered Newton equations we can apply all the methods and tools used 
for the homogeneous potentials that are described in Section |5l In particular, we have the 
following. 

Theorem 6.2. Assume that force F S (C;t-i[^])" has exactly D{n,k) proper Darboux points 
[d] G ©(f). Then A{d) satisfy the following relations: 

or, alternatively 

for < r < n — 1. 

This theorem has the same consequences for the Jacobi integrability of homogeneous 
forces as Theorem |5.4| for the integrability of homogeneous potentials. Namely, for a given 
k, there is only a finite set of candidates for A(d) satisfying the necessary conditions for 
the Jacobi integrability given by Theorem 16.21 Hence, we can try to perform an analysis 
similar to that for homogeneous potentials, and find all integrable forces for small k. Such 
an analysis was performed in ||4l|]. 

It appears that for n = 2 and k = 3 almost all forces F G C2[^] x C2[(j] are potential. 
For k = A and k = 5 several families of integrable forces were found. Among them some 
intriguing non-trivial examples appear, e.g., force with components 

Fi = qiqi, F2 = -^qiqi- 

is integrable in the Jacobi sense, with two first integrals which are both of degree four with 
respect to the velocities. They have the following form 

h = 24pip| + ?>qlm\vl + llqiqipipi - 3qlvi) + ^(>q\ql, 
h = 162pl{qip2 - qiVi) - ^li^lPi - IQqiqiVm + l^^lPi) + ^(^qWi- 
It is also worth to mention a remarkable family of forces given by 

Fi = \qiq\-^ F2 = q\-\ (6.12) 

where > 2 is and A G C. The Newton equations with this force admit the following first 
integral 

h = \v\'r\q\. 

If (fc. A) belongs to an item of the Morales-Ramis table l|4.10|l , then the system is integrable 
in the Jacobi sense with an additional polynomial first integral l2- Moreover, for an arbi- 
trary M > we find A such that the degree of I2 with respect to the momenta is greater 
than M, and there is no an additional polynomial first integral independent with l\ and 
having degree with respect to the momenta smaller or equal to M. Additionally, if (fc. A) 
belongs to an item different from 2 in table (|4.10|) . then there exist two additional poly- 
nomial first integrals I2 and I3 which are functionally independent together with Ji. The 



above statements were formulated in a form of a well justified conjecture in |4l[1 
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7 Open problems and perspectives 



Let us consider a two degrees of freedom Hamiltonian system with homogeneous polyno- 
mial potentials of degree k > 2. In ||39l] it was shown that if a potential of degree fc > 4 
admits a polynomial first integral of degree higher than two with respect to the momenta, 
then this first integral is a product of polynomial first integrals of lower degrees. One 
would like to prove the following. 

Conjecture 7.1. If a homogeneous polynomial potential of degree > 4 admits an additional 
polynomial first integral, then it admits an additional first integral of degree at most two with 
respect to the momenta. 

Let us consider a restricted version of this conjecture. Namely, let us assume that 
we consider only a generic potential of degree k > 5. Such potential has exactly k proper 
Darboux points, for each k we know two elements of its distinguished spectrum 'i^ ^, namely 

Ai^k={-l,-\,k-2,...,k-2), A2,k= (-^,k + l,...,k + l). 

^ V ' V 2k ^ ^ '/ 

k—2 times k—1 times 

Let us assume that for an arbitrary A: > 5 set 32 k only these two elements. Then 
we can prove the restricted version of the conjecture showing that the only potential with 
spectrum yi^- is the potential separable in the Cartesian coordinates, and the only potential 
with spectrum A2^k is a potential separable in parabolic coordinates. Unfortunately, the 
assumption about the number of elements of the set J2 k only for k < 13. For k = 14, 

set J2 k contains two additional elements 

A3,k = ( - 15, 15, 780, 5655) , 

7 times 

Au = ( - }Z_^, ^,^,39, 39, 5655) . 

8 times 

Then, additional elements of 3 2 k ^ppe^r for k = 17, 19,26,32, . . ., and it seems that there is 
no upper bound for these exceptional values of k. 

The presented method can be applied effectively only for small values of n and k. With 
the known computational algorithms for determination of distinguished spectra and recon- 
struction of potential, it seems that limiting values are n = 3 and k < 5. 

We believe that a substantial progress in this field will be possible after development 
of the higher order variational equations techniques. Till now, they are used only in cases 
when the first order variational equations have a very specific form, see js^. 
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